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Abstract 


We discuss general on-shell couplings of a scalar with two Z bosons using an operator analysis. 
In addition to the operator originated from the Higgs mechanism, two dimension-five operators, 
one CP-even and one CP-odd, are generated only at the loop-level. Simple formulas are derived 
for the differential decay distributions when the Z pair subsequently decay into four leptons by 
computing the helicity amplitudes, from which it is shown the CP-odd operator merely induces a 
phase shift in the azimuthal angular distribution between the two decay planes of the Z bosons. We 
also investigate new physics scenarios giving rise to loop-induced decays of a scalar into ZZ pair, 
and argue that the total decay width of such a scalar would be order-of-magnitude smaller than 
that of a Higgs boson, should such decays be observed in the early running of the LHC. Therefore, 
the total decay width alone is a strong indicator of the Higgs nature, or the lack thereof, of a scalar 
resonance in ZZ final states. In addition, we study the possibility of using the azimuthal angular 


distribution to disentangle effects among all three operators. 


I. INTRODUCTION 


In the standard model (SM), the electroweak gauge bosons obtain their masses through 
the Higgs mechanism, which postulates the existence of a scalar particle whose vacuum 
expectation value (VEV) breaks the electroweak SU (2); x U(1)y symmetry down to U(1)em. 
If the scalar, the Higgs boson, is a SU(2),;, doublet denoted by H = (h*,h)*, then its kinetic 


term 
2 


1 
s (1) 
contains mass terms for electroweak gauge bosons after the neutral component of the Higgs 


doublet gets a VEV, (H) = (0,v)"/./2, where g and g’ are the gauge couplings for the 


i ae 
Dir = (2, — ig W, — ig 5B, H 


SU(2)z and U(1)y, respectively, and o% are the Pauli matrices. Using the mass eigenbasis 
_ Wi sB, g, Wi + 9B, 


= 1 


v2 


one finds from Eq. (I) the following masses 


1 1 
my = 59, mz =5V9 +9", ma =Q. (3) 
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Furthermore, there are also three-point and four-point couplings from the Higgs kinetic term 


derived by replacing my > my(1 + h/v) in the gauge boson mass term: 


h 2 
(+4) mV, V”, (4) 


where V = W, Z. The form of the hV V coupling is completely determined by the electroweak 
gauge invariance to be 
E Ü (5) 
v 
Therefore, measurements of the three-point vertex in Eq. (5) will be a striking confirmation 
of the Higgs mechanism. 

Experimentally, the hVV vertex plays an important role in discovering the Higgs boson 
at the Large Hadron Collider (LHC). For a Higgs mass above 150 GeV or so, the branching 
ratio is dominated by decays into WW and ZZ ih In particular, h > ZZ — 42 is the gold- 
plated mode for the discovery of a moderately heavy (= 180 GeV) Higgs boson, which is a 
very clean signature with relatively small backgrounds. The excellent energy resolution of 


the reconstructed electrons and muons leads to a clear 4-lepton invariant mass peak, which 


allows for precise measurements of the mass and width of the Higgs boson |2]. 
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Given that so far all data from collider experiments agree with predictions of the SM 
quite well, there are very few experimental hints on what could (and could not) be seen at 
the LHC. Therefore, if a new scalar resonance is observed in the WW and ZZ final states, 
it is perhaps prudent to proceed without presuming the discovery of a Higgs boson whose 
VEV gives masses to the W and Z bosons. Only until after one could verify the decay 
indeed occurs through the three-point coupling in Eq. (5), can one gain some confidence in 
the Higgs mechanism as the origin of electroweak symmetry breaking (EWSB). 

In this work we study the physics giving rise to decays of a scalar into two Z bosons, 
with an emphasis on probing the Higgs nature of the scalar. Such a final state is interesting 
in its own right because of the high degree of symmetry in two identical spin-1 particles. 
Early studies of such systems resulted in the Landau-Yang theorem, which forbids decays 
of a spin-1 particle into two photons l. Recently similar arguments to the Landau-Yang 
theorem have been extended to decays of a massive spin-1 particle, the Z’ boson, into two Z 
bosons |4. There it was discovered that the azimuthal angle between the two decay planes 
of the Z is a very useful observable in discerning different interactions of the Z’ with the Z 
bosons. 

Here we consider the production of a scalar S in the gluon fusion channel, which is the 
dominant production mechanism of a Higgs boson at the LHC [1], and its subsequent decays 
into two Z bosons. We do not assume the scalar © plays the role of the Higgs boson in the 
Higgs mechanism. In particular, we point out non-Higgs-like couplings are induced only at 
the loop-level, and investigate in detail implications on the underlying new physics. Since 
we presume the scalar S and the Z bosons are all produced on-shell, implying mg > 2mz, 
our analysis is different and complimentary to studies on anomalous Higgs couplings in the 
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.) Differential 


vector boson fusion production, where the vector bosons are off-shell pi 
of anomalous Higgs couplings at the linear collider have been studied in 
distributions of a scalar decaying into ZZ — 4@ in the general case have been computed in 
Refs. |8|, (9. However, applying the symmetry argument as in Refs. l d would allow us 
to simplify the decay distributions dramatically, making it transparent the usefulness of the 


aforementioned azimuthal angle. We also argue that the total width of a scalar decaying to 


1 Tt is worth pointing out that at the LHC the production rate in the gluon fusion channel is an order of 
magnitude larger than the vector fusion production through out a wide range of Higgs mass. 


ZZ through loop-induced effects should be much smaller than that of a Higgs-like scalar, 
if the loop-induced decays should be observed at the LHC in the early running. Therefore 
measurements on the total width alone is a smoking gun signal for the Higgs nature of the 
scalar resonance. 

This paper is organized as follows: in the next section we compute the differential distri- 
bution of the decay of a scalar into ZZ — 4£ using the helicity amplitudes method, followed 
by a discussion on the possible new physics giving rise to loop-induced couplings. In Section 
IV we perform simulations on the total decay width measurements as well as azimuthal 
angular distributions between the two decay planes of the Z boson. Then we conclude in 
Section V. We also provide two appendices, one on a toy model in which the loop-induced 
coupling is mediated by the heavy W’-boson loop and the other on the Lorentz-invariant 


construction of the aforementioned azimuthal angle. 


II. HELICITY AMPLITUDES FOR S > Z(\1,k1)Z(A2, k2) > (€1£1) (262) 


We use the notation (A;, ki), i = 1,2, to denote the helicity state and momentum of the 
two Z bosons in the laboratory frame. Assuming all three particles are on-shell, the possible 
helicity states U2 of the Z pair are determined by conservation of angular momentum 
to be Ut*, U-~, and VW, from which we see the parity-even combinations are Y++ + Y~- 
and °° while the parity-odd one is Y+*+ — ~-~. In terms of effective Lagrangian, the three 


helicity amplitudes are described by the following three operators 


1 5 le " 1c V r7 po 
Leff = gms 5 (az Zy + mA Zw + Imz nena Zz ) ’ (6) 


where Z, = OZ, — 0, Z, is the field strength, and c;,i = 1, 2,3, are dimensionless constants. 
A fourth operator, Z,,,(Z"0"S — Z”O"S), is related to the cı and cz terms upon the equation 
of motion. The tensor structure of the decay amplitude of S + Z,(k*) + Zo(k8) is 
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ees Mag = Mg feh fo Jap — —> [Gap ki ` k2) — (k1)a(k2)g] + S cask te} s (7) 
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where ef and £ are the polarization tensors of Z, and Z3, respectively. Terms in Eq. (Ø 
proportional to cp and c3 are the so-called anomalous Higgs couplings. 


Following the method and convention in Ref. [4], we calculate the helicity amplitudes 


Mzz: 


2m? 4m? 
M - faa (1- T4) žic 1— a (8) 
Msg Msg 
m2, m 
Moo = ms fa (1 = 7s) + z . (9) 


Notice that Moo is real while the amplitudes M++ are complex in the presence non-zero c3. 


Therefore, we can parametrize the three helicity amplitudes in terms of two real numbers, 


Mr and Mz, and one phase 6: 
M4 = Mre”? 5 M= Mre~® ; Moo = ML, (10) 
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When the two Z bosons further decay into (¢;0;)(¢202), the phase 6 enters into the differential 
distribution in a simple way. To see this, recall that the angular distribution of the decay 
Z; > 44 in the rest frame of Z; has the dependence e*t: where m; = 0, +1 is the spin 
projection along the z axis and @; is the azimuthal angle. Since only the relative angle @ 
is physical we set ¢; = 0 and ġ = ¢. (See Fig. [I]}) From Eq. we see ô only enters as 
a phase shift in @é > ¢+ ô. Furthermore, the angular dependence of the differential decay 


rate is schematically 
— m |a + aet T® + age "Ot /? ~ by + bo coslo + 5) + bs cos(2 + 26). (13) 


For a Z’ boson decaying into the ZZ pair, the cos 2¢ term is absent in Eq. and a similar 
phase shift 6’ enters as ¢ > @+ 26’ i 


Using gz and gp to denote the coupling of the Z boson to the left-handed and right-handed 


leptons, respectively, we arrive at the differential distribution of S > Z,Z_ > (¢,6;) (los) 


FIG. 1: Two decay planes of Zı > 4&4 and Z2 > lək define the azimuthal angle ¢ € [0, 2x] which 
rotates l2 to lı in the transverse view. The polar angles 0, and 02 shown are defined in the rest 


frame of Zı and Z2, respectively. 
following the method of helicity amplitudes lid: 


dI 1 fl., 
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where the definition of 01,02, and @ are given in Fig. Integrating over the polar angles, 


we get the expression 


dv 1 [8 


T? Mr (92,- gN 16 / M? 
—— ô) + — | — +2 15 
+ > M (G8 cos(d + ô) + 9 (I2 (15) 
The normalization factor is given by integrating the above expression, 
32r / M? 
N = | — +2}. 16 
9 (3; t ) c 


Let’s consider turning on c; one at a time: 


ec, £0 and cp = c = 0: 


2 
ms 


=1- ae and 6=0. (17) 
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Mz 
Mr 


This is the case when S' plays the role of the Higgs boson in the Higgs mechanism. 


Since we assume Mms > 2mz for on-shell production, we see |Mr/Mr|az0o > 1 and 
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the longitudinal component of the Z dominates over the transverse components in the 


decay, especially in the limit of large mg. 
e co #0 and & = c3 = 0: 


Mz e E (5 


—1 
a (a and 6=0. (18) 
Mr co £0 1 = m2 /(2mz,) Mr D 


In this case |M;,/Mr|-.40 < 1 and the transverse polarization of the Z dominates in 


the decay. 


e c #0 and & = & = 0: 


My, 
Mr 


=U and ô = = ; (19) 
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This is a particularly simple case, as the normalized differential distribution in Eq. 


reduces to 
dV 1 1 
— =—(|l1-- 2 : 20 
Tdo =( pe 6) (20) 
Previous analysis assuming the SM Higgs boson can be found in Ref. hah, while Refs. 
addressed the CP violation due to the simultaneous presence of cı and c3 terms. Our 
general and simple result in Eqs. (14) and is in agreement with the lengthy expressions 


in Ref. [9]. Furthermore, our analysis makes it clear that the effect of a non-zero c3, which 


is CP-odd, is to induce a phase shift in the azimuthal angular distribution. 


Ill. NEW PHYSICS AND LOOP-INDUCED DECAYS OF S 


Among the three operators in Eq. (6), cı has the form of the three-point coupling in the 
Higgs mechanism and could be present at the tree-level,? while both c> and c3 are higher 
dimensional operators induced only at the loop level mi If cı = 0 at the tree-level, the 
scalar S is not responsible for giving W and Z bosons a mass. We also assume the existence 
of the following gluonic operators, 


Cg2 a a uv Cg3 vpo ya a 
P SG a” and TF per P G Gs (21) 


? cı could also be generated through dimension-five operators such as S|D,H|?, which is suppressed by a 


high mass scale comparing to Eq. (6). We will not consider this possibility further in this work. 


T 


so as to allow for the production of S in the gluon fusion channel. In the SM it is well-known 
that Cg2 is induced by the top-quark triangle loop when S is the Higgs boson. In fact, c2 
is also present in the SM through the W-boson as well as the top-quark loop [15], which 
is nonetheless overwhelmed by the tree-level cı given in Eq. (6). On the other hand, the 
CP-odd operators c3 and c,3 can be generated by a fermion triangle loop when the fermion 
has an axial coupling with the scalar S ; 

At the LHC, the event rate Bo(gg > S > ZZ) is 


(S > ZZ) 


r 
Bo(gg > S > ZZ) = o(gg > S) x Br(S > ZZ) =o(gg > S) x , (22) 


Ptotal 
where the total decay width Tota) is given by summing over all decay channels, including 
possible decays into SM fermion pair f f, 
Tro= XO I(S4VV)+S (S> ff). (23) 
V=g,W,Z,y7 f 


A few model-independent observations are in order: 


e While the decay channel into fermions may or may not exist, electroweak symmetry 
ensures the existence of decay channels into WW and yy once S —> ZZ is observed. 


Establishing the production gg — S also guarantees a decay channel into two gluons. 


e If the event rate Bo is comparable to the SM expectation of a Higgs boson, then the 


branching ratio into ZZ pair should be sizable 
Br(S > ZZ) > O(107"). (24) 


The SM Higgs production and decay gg > h > ZZ — 4 has an event rate in the 
order of 5 fb after multiplying o x Br with the pre-selection efficiency [2]. Therefore 
if ø x Br is an order-of-magnitude smaller that that of a SM Higgs, it would require 
an integrated luminosity of 300 fb~! to achieve 5ø significance for discovery, which is 


clearly beyond the early running of the LHC. 


e If S — ZZ is observed to occur through the loop-induced operators in the early LHC 
data, then a sizable Br(S — ZZ) implies the total decay width 


I(S > ZZ) 


eS 
ot Br(S > ZZ) 


(25) 


should also be one-loop suppressed. 


e Similarly, loop-induced S — ZZ and a sizable Br(S — ZZ) imply? 


riS 27) (Ss 


ca =| SOG). (26) 


2 
gi 

Since we expect c; and Cg; to be proportional to the electroweak coupling dey and the 
strong coupling a,, respectively, a large multiplicity factor (2 O(1)) in c; should be 


present. 


In the following we investigate new physics scenarios where the production and decay into ZZ 
of S occur predominantly through loop-induced operators. Such possibilities arise naturally 
if S is a SM singlet and couples to SM matter only through a messenger sector. In particular 
we focus on cases with a sizable branching ratio Br(S + ZZ) as in Eq. (24), so that S would 


have a comparable event rate to that of a SM Higgs boson. 


A. Fermion Loop-induced S —> gg 


In the SM gluon fusion production is induced by the top quark loop i, 


(SM) _ V2a, ms 
BO ae ag (27) 


It is well-known that this coefficient is related to the top contribution to the gluon two-point 
function from the Higgs low-energy theorem hd, If the messenger sector contains a pair 
of heavy vector-like fermions (Q°, Q) in the fundamental representation of SU (3)e with the 


interaction 
MQ + yaS QQ, (28) 
then its contribution to the gluon two-point function is 
1 TE E E 
=a 1- 1672 OF log 2 Gee, (29) 


where pe = 2/3 is the contribution to the one-loop beta function of QCD from a Dirac 


fermion and Mg(S) = mg + y@S is the mass of the new heavy fermion Q when turning 


3 Tn Eq. we have neglected an extra factor for decaying into massive gauge bosons, which is order unity 
unless mg is very close to the 2mz threshold. 


on the scalar as a background field S — S. To obtain the scalar-gluon-gluon coupling, the 
Higgs low-energy theorem instructs us to expand Eq. to the first order in S pi 


Cg2 = Z 7T YQ. (30) 


Strictly speaking, the low-energy theorem applies only when the mass of the particle in the 
loop is much larger than the scalar mass, m4 / (4mo) < 1, so that the loop diagram can be 
approximated by a dimension-five operator. We will always work in this limit in the present 


study. The partial width of S —> gg can be computed: 


3 
SS yj. (31) 


B. Fermion Loop-induced S > ZZ 


Next we consider the case when the messenger sector contains vector-like fermions (L°, L) 


charged under the electroweak gauge group with the interaction 
MLL L + y,S L°L, (32) 


where (L°, L) are in the fundamental representation of SU(2); and carry the hypercharge 
Yz under U(1)y. The contribution of L to the two-point function of the Z boson is simply 


1 2N, M? (S 
-2 i Pa (choy? + sib YRA) kg | Zua", (33) 


4| 16r s2 Y ie 
where po = 2/3, pi) = 4/3, Ne is the dimensionality of the SU(3),. representation L belongs 


to, and d® = 2 is the dimensionality of the SU (2) fundamental representation. In addition, 


Cw and s, are the cosine and sine of the Weinberg angle. Then we compute 


L Qem Ne (ct 42 A ms 
= 2s Yfd — 34 
Cg pym C2 82 Cw + Swt LOF m ( ) 
2 
1 
rO(S > ZZ) =P (=) ain (cy)’mg (35) 
S 
m2 a? N2 x 2 m3 
= P Z em c (c 264 Y 2q ’) ES. 2 36 
(=) 576T3 c4 s4, Cw + Sut DOr m2”: vo) 


where P(x) = vI — 4x (1 — 4x + 62) is a factor correcting for the massive final states in 
the decay width. Notice the additional difference between Eq. and Eq. due toa 


color factor of 8 since there are eight gluons in the final states for S —> gg. 
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It is worth commenting that, since SU (2)z is broken and gauge invariance does not forbid 
a mass term for the Z boson, one might expect a contribution to cı be generated at the one- 
loop level. However, recall that vector-like fermions do not give one-loop corrections to the 
Z boson mass term, and thus make no contributions to the m}Z,,Z" operator which would 
have given a contribution to cı after applying the Higgs low-energy theorem. This argument 


suggests that any contribution to cı at one-loop level would come from applying the Higgs 


low-energy theorem to operators with four-derivatives such as (O Z,,)*, which upon using 
the equation of motion is suppressed by (mz/mz)* and can be safely neglected for a heavy 
mz. We explicitly computed the fermion triangle loop diagram in a large mass expansion, 


my — oo, and verified that the first contribution to cı indeed starts at (mz/my)*. 


C. Gauge Boson Loop-induced S > ZZ 


The last possibility we consider is when the messenger sector includes a new set of heavy 
electroweak gauge bosons (W’, Z’). In the SM the W contribution to the loop-induced 
decay of the Higgs into yy dominates over the one from the top-quark loop due to a large 
beta function coefficient “7” multiplying the W loop result [16]. One may expect a similar 
situation for the W” contribution to the loop-induced decay into ZZ. Given the existence of 
two sets of electroweak gauge bosons, the simplest model must contain two copies of gauged 


SU(2). Schematically, the symmetry breaking pattern is a two-step process: 
SU(2), x SU(2)2 x U(1)y > SU(2), x UA)y > U (1)em , (37) 


where the two SU(2)’s are broken down to the vectorial subgroup, identified with SU(2)_, 
at a high scale fı using a linear sigma model. Subsequently SU(2); x U(1)y is broken down 
to U(1)em at a low scale fọ = v following the Higgsless model id 

In the Appendix A we explicitly construct a toy model whose gauge sector is the same 
as the three-site Higgsless model [dbi although we are interested in a different corner of 
parameter space, € = (f2/ f1)? < 1. For example, if fı ~ 1 TeV and fọ = v & 246 GeV, 
we have e ~ 0.06 < 1. In this case the W’ and Z’ can be as light as several hundreds GeV 
for weakly coupled theories. We computed the W” contribution to the ZZ self-energy. At 


leading order in €, 


e? MÈ, (S) ° 
— —— (7ci,) log a Lge, (38) 


where we see the same large coefficient as in the scalar coupling to two photons. In the set 


up we have in the Appendix A, Mj,,(S) = (9? + 93)(f1 + S}? at leading order in e, which 


leads to 
+ em M 
W em Cw S 

= a 39 
C2 On s2, fi ( ) 

Using Eq. (35) we arrive at 

i 2\ 4902 c4 m3 

Preuss ZZ)=P mz em w S 40 
( ) te | 250n° s4, f? ey) 


There is also a contribution to cı induced at one-loop by the W’ boson that is suppressed 


by e, which we ignore. 


Given Eqs. BI), (85), and (40), we can now compare T(S —> gg) with T(S > ZZ) and 
see that the decay width into ZZ can easily be comparable to the decay width into two 
gluons. This is especially the case for the W” loop due to the large coefficient in Eq. (89): 


(w’) 2 2 
mG > ZZ) VIP (=) BACA (41) 


r(S —> gg) $ 


Even in the case of purely fermionic contribution in S + ZZ, assuming the fermion (L°, L) 


carries no hypercharge, the ratio of the two partial widths is 


(f) 2 2,2 
rY > ZZ) 0.01N2 x P (z) ee (42) 
TS > gg) mg] MYO 


which could still be O(0.1) if the multiplicity factor Ne 2 3. This could be achieved if the 
fermion (L°, L) is also in the fundamental representation of SU (3)e, resulting in N. = 3. In 
the end, we have demonstrated that new physics scenarios giving rise to loop-induced decays 


of S into gauge bosons could easily give a significant branching ratio into ZZ bosons. 


IV. OBSERVABLES AND SIMULATIONS 


In this section we discuss two observables which could be useful in disentangling whether 
the scalar coupling to Z bosons is as predicted by the Higgs mechanism or induced by new 


physics at the loop-level. 
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FIG. 2: The dashed line is the total decay width for a SM Higgs boson and the solid line is that 
of a scalar S whose width is three orders of magnitude smaller. The yellow (shaded) region is the 


detector resolution. 


A. The Line Shape 


Among the three possible on-shell couplings of S with Z bosons, only cı could be present 
at the tree-level with an order unity coupling when S plays the role of the Higgs boson in the 
Higgs mechanism, while both c2 and c3 are non-zero only at one-loop level. This observation 
suggests that the total width of a scalar decaying through cy and c3 must be much smaller 
than that of a scalar decaying through cı, in order for the decay channel to be observable 
in the early LHC running. Using the SM Higgs as an example, the partial decay width 
T(h > VV) is 


m% \ Grm? 
yP | — ua V=W,Z, 43 
E (Fr) 16/27 (28) 


where Gp is the Fermi constant and ôw = 267 = 2. Comparing with the partial decay width 
of a W’-loop induced decay, we see 


TW) (S => ZZ) 


~ 1072 44 
Tih ZZ) uae (44) 


for a W’ mass as light as 500 GeV. As emphasized previously, in order for the event gg > 
S — ZZ to be observable at the LHC with say 30 fb~! luminosity, the branching ratio 
Br(S — ZZ) should be sizable and comparable to that of a SM Higgs into ZZ. It then 
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FIG. 3: The ZZ invariant mass distribution for a SM Higgs boson and a scalar S decaying through 
loop-induced effects, using a 2 GeV bin size. The narrow width of S is below the detector resolution, 
resulting in a concentration of all events in just one bin. Note that for a sufficiently small bin size 
one would resolve the peak, albeit with a form which is dominated by the detector resolution (a 
Gaussian, if the usual assumptions of detector smearing are made). In the plot we assume the 


event rate of gg > S —> ZZ — Al is only 10% of rate for the SM Higgs boson. 


follows that 
Tot (S) 7 r(S > ZZ) Br(h > ZZ) 


Math) BOs) TSZ) 


In other words, we would observe an extremely narrow peak in the ZZ invariant mass 


~ 10°. (45) 


spectrum if the scalar S only decays via cy and c3. In fact, the peak is so narrow that the 
width is completely below the detector resolution. In this study we use a 2 GeV bin size 
which is comparable to the energy resolution of the detector. On the other hand, a scalar 
participating in the EWSB like the Higgs boson would have a width above the detector 
resolution at the LHC except when its mass is below 200 GeV. Therefore, the Breit-Wigner 
line shape of a scalar resonance in the ZZ invariant mass spectrum is a strong indicator on 
the Higgs nature (or the lack thereof) of the scalar. 

In Fig. P] we show the total decay widths of a SM Higgs boson and a scalar S decaying 


through loop-induced operators, and compare with the detector resolution at the LHC using 
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the following lepton energy smearing: 


ôE a 


2 —— © b , 

E ./E/GeV 
where a = 13.4%, b = 2%, and ® denotes a sum in quadrature 29}, We see while the 
width of the SM Higgs could be resolved above a 200 GeV mass, the small width of the 


(46) 


S is completely buried in the detector resolution. The narrow width of S implies, in the 
invariant mass distribution of the two Z bosons, all the events would be concentrated in 
just one bin, resulting in a spectacular resonance peak even if the event rate is smaller than 
that of a SM Higgs. However, for a sufficiently small bin size, one would resolve the peak, 
albeit with a form which is dominated by the detector resolution (a Gaussian, if the usual 
assumptions of detector smearing are made). In Fig. [3] we simulate the ZZ invariant mass 
distribution for a SM Higgs and the S scalar using a 2 GeV bin size. To be conservative, in 
the plot we assume the event rate Bo(gg > S > ZZ — 4€) is only 10% of the SM Higgs. 
It is then clear that the total width measurement would allow for a distinction between the 
Higgs and a scalar © which decays only at the loop-level, except when the Higgs has a mass 


below 200 GeV and its width is comparable to the detector resolution. 


B. Angular Distributions in ¢ 


In the following we consider the dependence on the azimuthal angle between the two 
Z decay planes in the normalized differential rate in Eq. (15), by turning on one operator 
at a time. By considering the normalized rate, the dependence on the magnitude of the 
coefficients c; in Eq. (6) drops out and the angular dependence is largely determined by 
kinematics. It is worth mentioning that in Eq. the cos(@+ 6) term is highly suppressed 
due to the approximate symmetry g? ~ gẹ in the leptonic decays, so only the cos(2¢ + 26) 
term and the constant term will contribute. As can be seen from Eq. (17), for a Higgs-like 
scalar, cı # 0, the constant term becomes more dominant as the mass gets larger. On the 
other hand, for cp 4 0 the cos 2¢ terms is more important for a heavy scalar. 

In Fig. Ø| we simulate the azimuthal angular dependence in the normalized decay distri- 
bution for two different scalar masses: 200 and 400 GeV. To facilitate future experimental 
analysis, we provide in Appendix B Lorentz-invariant expressions for various angular vari- 


ables in the decay into ZZ — 4@, including the azimuthal angle ¢. We select the events by 
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FIG. 4: The normalized azimuthal angular distributions for 200 and 400 GeV scalar masses, turning 


on one operator at a time. 
requiring the following cuts on the lepton transverse momentum pr and rapidity 7: 
pr > 15 GeV and In| < 2.4. (47) 


In the simulation we assume the SM background coming from qq > ZZ — 4¢ has been 
reduced in the data sample using existing procedures for searching for a SM Higgs boson [|2]. 
The angular dependence of the SM background and its interplay in the Higgs search was 
studied in Ref. [23]. In the case of loop-induced decays, one can take advantage of the 
extremely narrow width and impose stringent cuts on the ZZ invariant mass to reduce the 
backgrounds. Therefore we do not include backgrounds in the plot. 

From Fig. [4] we see that the CP-odd case can be distinguished from the CP-even case 
in the angular distribution, which has been discussed in Refs. {11}, |24. The comparison of 
tree-level versus loop-induced operators in the CP-even case, however, does not seem to exist 
in the literature, to the best of our knowledge. We see that, even though cı 4 0 and cp Æ 0 
have the same phase in the angular distribution, the magnitudes are different even for a 
low mass of 200 GeV. Recall that this is also the mass range where the width measurement 
could be biased by the detector resolution. So one could use the angular distribution as an 


extra handle to distinguish a Higgs boson from a non-Higgs-like scalar. 
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V. CONCLUSION 


In this work we considered the most general on-shell couplings of a scalar with two Z 
bosons. In the SM the decay of the Higgs boson into ZZ final states is the gold-plated mode 
for discovery due to excellent energy resolution for charged leptons. However, in order to 
verify the Higgs mechanism as the origin of mass for the electroweak gauge bosons, it is 
necessary to measure the coupling between the scalar and the gauge bosons. By using an 
operator analysis for the most general couplings, we point out that dimension-five operators 
responsible for the anomalous Higgs couplings are generated only at the loop-level, while 
the Higgs mechanism would lead to a dimension-three operator at the tree-level. 

Using the method of helicity amplitudes, we computed the differential decay distribution 
of a scalar decaying into ZZ — 4. Our formulas are simpler than and agree with previous 
calculations. Furthermore, our results make clear the advantage of using the azimuthal 
angle between the two decay planes of the Z bosons in discerning effects between CP-odd 
and CP-even operators. 

If the scalar is produced in the gluon fusion channel, which gives the largest production 
cross-section for the Higgs boson at the LHC, a decay channel into two gluons must also 
exist. Then in order for the event gg > S — ZZ to be observable at the LHC in the 
early running, the partial decay width should be comparable to the partial width into two 
gluons. We investigated new physics scenarios giving rise to such a possibility by considering 
fermion-loop and W’-loop induced couplings of a scalar with ZZ bosons. 

One important implication of loop-induced operators is that the total width of a non- 
Higgs-like scalar, if its decay were discovered at the LHC early on, should be order-of- 
magnitude smaller than that of a Higgs-like scalar, which decays through tree-level processes. 
Again this is a corollary of requiring a sizable branching ratio into ZZ final states from loop- 
induced effects. Therefore measurements of the total width of a scalar resonance in final 
states with two Z bosons is a strong indicator on the Higgs nature of the resonance, except 
when the scalar mass is below 200 GeV and the SM Higgs width is comparable to detector 
resolution. In this regard, azimuthal angular distribution could provide an extra handle 
in determining not only the CP property of the scalar but also whether the decay is loop- 
induced. Only when the scalar coupling with the Z bosons is verified to be the one as 


predicted by the Higgs mechanism, can one gain confidence in the Higgs mechanism as the 
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origin of electroweak symmetry breaking as well as the discovery of a Higgs boson. 
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Appendix A: ZZ Self Energy From the Heavy Gauge Boson Loop 


In this Appendix, we compute the one-loop corrections to the Z self energies which are 
needed to construct the SZZ effective coupling using the low-energy Higgs theorem. For 
concreteness, we consider a simple gauge extension of the Standard Model (SM) which is 
based on the gauge group SU(2),; x SU(2)2 x U(1)y. We use two link fields X; and X», 

S a O_O 
Denn. (=i (48) 


1 
Eo = e20/f , (49) 


uy 


which transform as bidoublets under SU(2), x SU(2)2 and SU(2),; x U(1)y, respectively. 
The remaining unbroken gauge group is identified with U(1).m, whose generator is Q = 
T. (t) + Te + Y/2. Notice that the gauge sector of this model is identical to the so-called 
three-site Higgsless model studied in Refs. hkdh, except that we are allowing for a scalar 
degree of freedom in the radial excitation of X4. More importantly, we are interested in the 
limit € = (f2/ f1)? <1, which is also different from the three-site Higgsless model. 


The covariant derivatives are written as 


2 0" a . o° a 

DX = ô% — igy Wir + go> Woy (50) 
NN o ae ef" 

D,X2 = X2 — igaz Woy%2 + ideg z Pm (51) 
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where Wf, and g; are the gauge fields and coupling strengths belonging to SU(2);, i = 1, 2, 
respectively. Similarly B,, and g’ correspond to the gauge field and coupling of the U(1)y. 
The gauge bosons in the model obtain masses through the kinetic terms 
fF (p,34)"( DYE Bop yy)" DS 52 
oa p%1)'( e 22)" 2), (52) 


which lead to the the mass matrix for neutral gauge bosons (W?, W3, B): 


gf? -gg ff 0 
3 —ngof? B+ f2) dR (53) 
0 —g9'g f3 PF 


This matrix can be diagonalized by means of an orthogonal matrix which we shall call R: 


wÈ, Aj 
Wey = Ri Zu ? (54) 
B, Z, 

where the mass eigenstates are denoted by A, Z, and Z’. The eigenstate A is massless and 


identified as the photon. The couplings of our theory are related to the electric charge by 
- (55) 


cos Oy 


e 1 


e 
qh = oh Q= SSS 
cos ġ sin Oy ” sin dsin Ow ’ 


where Oy is the weak mixing angle (in the limit € > 0) and ¢ is an additional mixing angle. 


The other two eigenmasses are 

my = fe +9) ji- (56) 

2 (gi + 92) 
4 

my = fal D 1- i > (57) 

2 (gi + 93) 
where we have dropped O(e”) terms and 
1 1 1 


Clearly, Z is identified with the SM Z boson while Z’ is referred to as the heavy Z boson. 


cos ġ sin Ow sin ¢ sin 0w cos Oy 
3 «2 : 4 
R = | cos¢cosOw + e€ sos e a 2 sin dcos Ow — e a Z — sin Ow (59) 


For small e, the mixing matrix R has the following approximate form: 


— sin ¢ + esin ọ cost d cos¢+esin?¢cos?d —e tan Ow sin ġ cos? d 
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from which it is simple to verify that tree-level couplings between the scalar S and the Z 
boson start only at order e?. In other words, in this model cı = 0 at tree-level if we only 
keep terms up to O(e). The charged gauge boson sector can be worked out in a similar way, 
where the light mass eigenstate is identified with the SM W boson and the heavy eigenstate 
is denoted by W’. The couplings between S and the W’ and Z’ bosons have the form as 
predicted by the Higgs mechanism: 


S 2 
mi, (1 + 5) Vive, V=W,Z , (60) 


which is valid at leading order in e€. 

We would like to compute the one-loop correction to the Z self energy arising from 
loops of W’ gauge bosons. Unfortunately, much like the analogous corrections in the SM, 
the corrections to the two-point functions depend non-trivially on the particular Re gauge 
used to define the W” propagator |25]. However, by extracting Re gauge-dependent pieces 
from other one-loop corrections (i.e., vertex and box corrections) and summing these with 
those from the two-point function one can obtain an expression which is independent of the 
particular gauge chosen to do the calculation. This method, which is known as the Pinch 
Technique (PT), has been applied to the SM to obtain gauge-independent expressions for 
the gauge boson self-energies [26]. More recently, though, it has been extended to models 
with extended gauge sectors such as the model considered here . In this work, we 
will directly apply the results from the above references by taking the limit of our interest, 
e = (fa/ f)? & 1. We refer interested readers to Ref. for details and only make 
the following two comments. First, our results are obtained by taking the so-called “ideal 
localization” limit for the de-localized fermion introduced in Ref. “ Such a limit has the 
advantage of reducing the tree-level S parameter in the model. However, the main reason 
in our case is to decouple the de-localized fermion from the W’ boson, so as to remove the 
extra pinch contribution to the two-point function that is unnecessary for maintaining the 
gauge invariance. Second, even though we are allowing for a scalar degree of freedom in the 
radial excitation of X1, which is absent in the three-site Higgsless model, the computation 


in Ref. still carries through because S has no couplings to the Z boson in the order we 


4 Tn our model ideal localization is achieved by choosing the delocalization parameter x1, which is defined 
in Ref. leo, to be sin? ¢/(sin? ¢ — cos? ¢). 
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are working. Therefore the Z self energy in the non-linear sigma model (Higgsless model) 
is the same as in the (partially) linear sigma model we consider. 
The one-loop expression for the Z self-energy computed using the PT are then given by 


(up to O(e?)): 


w! Qem 3 amy f 
I’ (p?) Ee a cos® ġ sin? o 


2 


A 
+p (T, — 14€ cos? dcos2¢ c) log ——, 
My 


(61) 


where A is the cutoff of our effective theory. Notice that formally (mw /mw)? ~ 1/e so 
the longitudinal piece is considered mj,, x O(c), while the leading term in the transverse 
component has a large coefficient “7,” the same as in the SM W contribution to the photon 


self energy, which is to be expected. 


Appendix B: A Lorentz-Invariant Construction of ¢ 


In this Appendix we provide a Lorentz-invariant expression for the azimuthal angle @ 
between the two decay planes of the ZZ pair, so as to facilitate the analysis of angular 
distribution in ¢. Let pı and p be the momenta of the lepton pair coming from one Z, and 
p3 and p4 be the momenta of the lepton pair from the other Z. The parent momentum is 
P = pı + p2 + p3 + pa, which satisfies the on-shell dispersion relation P? = M?. We follow 
the notation in Ref. i pı = £1, po = h, p3 = 2, pa = b2. (See also Fig. M) 

In the rest frame of P, our azimuthal angle ¢ is given by 


Pi X P2 P3 X Pa 


m . ————— = cO, (62) 
[Pi|[P2|sin 12 |ps||p4] sin 34 


where the triple products in the numerator can be written in a Lorentz-invariant fashion: 
(Pi x pa) = Se” pip, P, = er, (pg x pai = LeriP, (63) 
M ee M M 
Note that we define e!?30 = 1 = egi93 = —e°!*3. Then it follows 


. > 5 P 
pıpouP -pspavP Pı: P3 Pı:Pa Pı 
€ E€ 


P -p3 P- p, M? 
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To arrive at a covariant expression for Eq. (62), we need to cast the denominator in the 
covariant form as well: 
2 


1 = 
Ipil|=—m-P, cos #49 = 1 mi 


— ———_ , 65 
M 2|pı||p2] l ) 


where m7, = (p; + p;)*, and similarly for |p2|,|ps|, and |p4|. In the end we have 


Pir P3 Pepa Dre LP 
M? | p3 -p3 p2- pa p2-P 
P -p3 P- p, M? 
cos g = = A r] 
M2m? M2m2 
(Pi P)(pa- P)(ps- P)(pa P)af1 — (1 E aes | an (1 = | 
(66) 
On the other hand, sin ¢ can be evaluated by the following relation, 


1 ¢P1p2pspa [Pi + Pol 


sin @ = = A 
M |p,| |p2| [ps] [pa] sin 612 sin 834 


(67) 


where 
Ey pt pi ri 
_ Ey p3 p3 D3 
= Europ Ph Py PS Di = aa = —Mp3 X pa: P1 - (68) 
Es P3 P3 P3 
E, pi pi Pi 


¢P1P2P3P4 


The covariant form is given by 


1 
: 1 M4) eP1P2P3Pp4 
sin ọ = = 


(Pi P)(pa : P)(ps - P)(pa - P)yjd — (1 = | (1 7 Sis] 


with A = 1 + mi,/M* + m3,/M* — 2m?./M? — 2m3,/M? — 2m?,m3,/M? . 
We can also determine the polar angle of pı in the rest frame of the 12 pair. A simple 


Lorentz boost gives 


_ E+E 
E = yE (1 + 8 cos 8) = ie N , (70) 
M192 2 
which leads to S 7 
Ey- E 2 
cos ĝ = ———" = (pı: P — p2: P). (71) 


[pi + Pol MA 
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For the polar angle of p3 in the rest frame of the 34 pair, simply replace pı and pə by p3 and 


pa, respectively, in the above. 
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